Hedin's equations and enumeration of Feynman diagrams 
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Hedin's equations are solved perturbatively in zero dimension to count Feynman graphs for self- 
energy, polarization, propagator, effective potential and vertex function in a many-body theory of 
fermions with two-body interaction. Counting numbers are also obtained in the GW approximation. 
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I. INTRODUCTION 

Consider a system of fermions described by the Hamil- 
tonian H = Y2i M*) + "52i<j v (hj)> where h is a single- 
particle operator and v is the two-particle interaction. 
In general, to obtain dressed correlators that describe 
quasi-particles and collective excitations, one must sum 
over classes of diagrams of all orders of perturbation the- 
ory in v. The random phase approximation for the ef- 
fective potential or the ladder expansion for the vertex 
and the T-matrix are two well known examples^. It- 
erative schemes are used in conserving approximations 
for the self-energjiSi^. The perturbative treatment of the 
electron gas requires a mixed approach, with partial re- 
summations; a late step in the long history of the r s - 
expansion of correlation energy is the evaluation of third 
order diagrams^. In these cases, it might be useful, or at 
least interesting, to know in advance how many diagrams 
are required at each order of the approximation, and how 
many are left out. The counting problem is much simpler 
than the original one, and is obtained by translating the 
recursive equations at hand to zero dimension. Proper- 
ties other than loop content or perturbative numbers may 
also be evaluated: a notable example is diagram enumer- 
ation according to their topologies, that can be mapped 
to models of statistical mechanics on random lattices^. 

In this report, I show how Feynman diagrams can be 
enumerated for the exact theory and its GW approxima- 
tion. A useful starting point is the set of five Hedin's 
equations 6 for the propagator (2(1,2), the effective po- 
tential W(l, 2), the proper self-energy S(l, 2), the polar- 
ization 11(1, 2), and the vertex T(l; 2, 3): 



G(1,2)= 5 (1,2)+ 5 (1,1')£(1',2')G(2',2) 
W(l, 2) = v(l, 2) + «(1, l')n(l', 2')W(2', 2) 
S(l, 2) = iT(2'; 1, l')G(l', 2)W(2', 2) 
n(l, 2) = -2ir(l; 2', 2")G(2, 2')G(2", 2) 
r(l; 2, 3) = 5(1, 2)tf(l, 3) + 

<5£(2,3) 



r(l;2',3')G(2",2')G(3',3"; 



<5G(2",3") 



(1) 
(2) 
(3) 
(4) 

(5) 



Repeated primed space-time and internal variables are 
hereafter understood to be summed or integrated. g(l, 2) 
is the Green function of the interacting system in the 
Hartree approximation, with exact particle density, so 



that Hartree-type insertion (tadpoles) are already ac- 
counted for. In the vertex equation, the Bethe-Salpeter 
kernel is specified as a functional derivative in the same 
quantities involved. Therefore, the five exact Hedin's 
equations are formally closed. A different functional clo- 
sure of Schwinger-Dyson equations, leading to diagram- 
matic recursion schemes, was developed by Kleinert et 

ai&a. 

It is useful to restate eq.(0 with the Hartree propaga- 
tor g replacing the exact one, G: 



T(l; 2, 3) = 6(1, 2)6(1, 3) + <?(!', !)<?(!, 1") 



<5S(2,3) 
5g(V,l") 



(6) 



It states that vertex diagrams are obtained from self- 
energy ones through the insertion of a bare vertex in a 
Hartree propagator line, in all possible ways. A deriva- 
tion is given in the end. The formula can be advantageous 
when the self-energy is approximated by some explicit ex- 
pression of g and v, which are the free functional variables 
of the many body problem. 

To cope with the vertex equation is a formidable 
t;is]«^ili. A strong non-trivial simplification is the GW 
approximation (GWA), where all corrections to the bare 
vertex r(°)(l; 2,3) = 6(1, 2)6(1, 3) are ignored. Eqs. <Q} 
0}, with To replacing T, then become a closed set of in- 
tegral equationsiSiiii. 

The upmost difficulty of the problem of solving Hedin's 
or GW equations, vanishes in zero dimension of space- 
time, where the perturbative solution in the interaction 
v is a mean to count Feynman diagrams. The con- 
straint of no Hartree-type insertions makes this approach 
much simpler than the usual one, based on the functional 
integraii^iiii. Counting numbers are also evaluated in the 
GWA, and the omission of vertex corrections makes the 
number of diagrams grow with a power law, instead of 
factorially. 



II. COUNTING FEYNMAN DIAGRAMS 

In zero dimension of space-time, the four Hedin's equa- 
tions (111 111 become algebraic, with variables g and v, and 
the functional derivative in the vertex equation JB] ) is 
an ordinary one. After removing imaginary factors and 
replacing the factor (—2), due to the fermionic loop and 
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spin summation, with a parameter I that counts the same 
loops, Hedin's equations are: 



G = g + gEG, W = v + vUW, 



(7) 



,0X3 



X3 = GWT, T1 = £G 2 T, T = l+g 2 — (8) 

By searching solutions as series expansions in the vari- 
ables v and £, one obtains coefficients that count all Feyn- 
man graphs, with weight one, that contribute to a per- 
turbative order (specified by the power of v), with a given 
number of fermionic loops (the power of £) . I begin by 
solving for the self-energy. Elimination of G, W and II 
gives 



X3(l - , 9 X3) 2 = gvT[l - (1 - £) 9 E] 



(9) 



It is useful to consider adimensional quantities, S = 
gvs(x) and x = g 2 v. After use of the equation for T 
one obtains a differential equation for s(x) 

ds 

a(l - is) 2 = (1 + xs + 2x 2 — )[1 - xs(l - £11 (10) 

ox 

There is no standard solution of this (Abel) equation. 
The perturbative expansion is evaluated 

E/vg = 1 + (2 + £)x + (10 + 9£ + £ 2 )x 2 + 
(74 + 91£ + 2U 2 + £ 3 )x 3 + 

(706 + 1063£ + 416£ 2 + 46£ 3 + £ 4 )x 4 + (11) 
(8162 + 14193£ + 7344£ 2 + 1350£ 3 + 80£ 4 + £ 5 )x 5 + .. 
= 1 + 3a; + 20x 2 + 189a; 3 + 2232.x 4 + 

31130a; 5 + .. (12) 

The last line corresponds to t = 1. For example, in eq. 
(|12fl we read that at order v 3 there are 20 self-energy 
diagrams; eq. I|ll|) gives the more detailed information 
that 10 graphs come with no fermion loop, 9 with a single 
loop and 1 with two loops. At each perturbative order 
there is a single diagram with largest number of loops; 
the sum of such diagrams yields the ring approximation 1 
S r = igW r , where the effective potential W r is evaluated 
with ring insertions II r = —2ig 2 . 

The expansions for the vertex and the polarization are 
then obtained (I omit the propagator and the effective 
potential) 

T = 1 + x + (6 + 3£)x 2 + (50 + 45£ + M 2 )x 3 + 

(518 + 637£ + 161£ 2 + 7£ 3 )x 4 + 

(6354 + 9567£ + 3744£ 2 + 414£ 3 + 9£ 4 )a; 5 + .(13) 
= 1 + x + 9x 2 + 100a; 3 + 1323a; 4 + 20088a; 5 + .(14) 

U/g 2 £ = 1 + 3a; + (15 + 5£)x 2 + (105 + 77£ + 7£ 2 )x 3 + 

(945 + 1044£ + 234£ 2 + 9£ 3 )x 4 + 

(10395 + 14784£ + 5390£ 2 + 550£ 3 + ll£ 4 )a; 5 + .. (15) 
= 1 + 3x + 20a; 2 + 189a; 3 + 2232a; 4 + 31130a; 5 + .. (16) 



Note that the counting numbers for the 1=1 expansions 
(|12|l and (|16fl of the self-energy and the polarization are 
the same. This occurs also in Q.E.D. though with smaller 
counting numbers, due to the cancellation of pairs of di- 
agrams which differ by orientation of a fermionic loop 
with an odd number of propagators (Furry's theorem) 16 . 

I now evaluate the counting numbers in the GWA, 
which corresponds to T = 1 in eq.©. The equation for 
the self-energy becomes algebraic cubic, with solution 

E GW /vg = 1 + (1 + t)x + (2 + U + £ 2 )x 2 + 
(5 + 15£ + 9£ 2 + £ 3 )x 3 + 

(14 + 56£ + 56£ 2 + 1U 3 + £ 4 )x A + (17) 
(42 + 210£ + 300£ 2 + 150£ 3 + 25£ 4 + £ 5 )a; 5 + .. 
= 1 + 2a; + 7a; 2 + 30a; 3 + 143a; 4 + 728a; 5 + .. (18) 

For £ = 1, the cubic equation is s(l—sx) 2 = 1. By solving 
for the inverse function, x(s) = s _1 — s~ 3 / 2 , one locates 
a singular point for s'(x) in x c = 4/27. The finite radius 
of convergence for the perturbative series i|18|) in GWA, 
implies that the number of self energy diagrams of order 
n grows with the power law (27/4)™. In GWA the vertex 
function is trivial, and the polarization Hgw — ^iG^Gw 
is the GW-dressed ring diagram: 

Hcw/g 2 £ = 1 + 2x + {5 + 2£)x 2 + 

(U+U£+2£ 2 )x 3 + 

(42 + 72£+27£ 2 + 2£ 3 )a; 4 + 

(132 + 330^ + 220£ 2 + U£ 3 + 2£ 4 )x 5 + .. (19) 
= 1 + 2a; + 7a; 2 + 30a; 3 + 143a; 4 + 728a; 5 + .. (20) 

The £ = 1 expansions for the self-energy and the polar- 
ization are again the same, as one can show that H/g 2 
and s solve the same cubic equation. 

Conclusion: By translating to zero dimension a closed 
system of equations for many-body correlators, one may 
enumerate the Feynman diagrams involved. I have shown 
this for the full theory of interacting fermions and the 
GW approximation, with the rule that a line is a Hartree 
propagator. 
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APPENDIX A 

I show that the vertex equation (jSJ) can be written in 
the form ©; the simple proof is a modification of Hedin's 
construction, to which I refer. I recall the definition of 
the vertex function 



r(l; 2, 3) = 6(1, 2)6(1, 3) 



(5X3(2,3) 
6V(1) 



(Al) 
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The potential V is the sum of an external potential 
(which is turned off at the end) and the Hartree poten- 
tial / d2v(l,2)n(2), built with the exact particle density. 
Hcdin's formula is obtained by writing the functional 
derivative in V, through the chain rule, as a derivative 
in the exact Green function G. It is admissible to do so 
in terms of g: 

5E(2,3) _ 511(2,3) 8g{V, I") 
SV{1) ~ Sg(V, 1") SV(l) 



From the Hartree equation (id(ti) — h(l') — 
V{l'))g{l', 1") = 5(1', 1"), one evaluates the required 
functional derivative 5g(l',l")/SV(l) = g(l', l)g(l, 1"). 
This ends the proof. 
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